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1 Introduction 

Consider an elastic thin membrane which occupies a planar region represented by a simply 
connected bounded open set ft C R 2 . In the same plane a vector force field T is applied on its 
boundary 9ft, so that the membrane is subject to a distributed pretraction state, expressed by a 
positive and symmetric tensor a = {cry }, i, j = 1,2, which satisfies the plane equilibrium equations 

diver = , in ft , (1.1) 
gv = T , on 9ft , (1.2) 

where v is the outer unit normal to 9ft. The transverse displacement u(x) of the membrane will be 
governed by the equation 

-div(o-Vu) = / , in ft , (1.3) 
u = ip , on 9ft , (1-4) 

where f{x) represents the distributed transverse load applied to it, and tp represents the prescribed 
transverse displacement at the boundary. 

In this note we wish to investigate the inverse problem of determining the plane traction state 
tensor a from boundary measurements on displacements and on the corresponding forces. As an 
initial attempt, since the external load / has no influence on the pretraction state a, and in analogy 
with many other well-known inverse boundary problems, see for instance [H[51[T5], it seems natural 
to treat the case when / = in (jl.3p . and consider as available data an arbitrary transverse 
displacement tp on the boundary and the corresponding transverse load on the boundary, namely 
the reaction of the boundary constraints 

A a tp = erVit ■ v . (1-5) 

Hence, fixing any tp £ _ff 1/,2 (9ft), if we denote by u £ _ff x (ft) the weak solution to the Dirichlet 
problem 



div(o-Vu) = , in ft , (1.6) 

u = tp , on 9ft , (1-7) 

we introduce the Dirichlet-to-Neumann (D-N) map as the bounded linear operator 

A a : Hi(dn) ^ H-i(dn) , (1.8) 

defined, in weak terms, by the formula 

< A (7 tp,V\ dn >= / crVu • Vw , for every v £ H 1 (ft) , (1.9) 
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where viqq denotes the trace on dfl of v G 

Thus we examine here the problem of determining a from the knowledge of A CT . The peculiarity 
of this problem is that, by its own nature, the tensor a is anisotropic, and since Tartar's example, as 
reported in [11) , it is well-known that a general anisotropic tensor a cannot be uniquely determined 
by the D-N map A CT . In this case however, we shall see, in the next Section [2] that we can take 
advantage of the null divergence condition (|l.ip and obtain the uniqueness in such restricted class 
of genuinely anisotropic tensors, see Theorem 12. 11 

In Section [3] we also consider the question of stability, that is of the continuous dependence 
upon the data. And although the present results are very preliminary, they show up interesting 
phenomena, which are markedly different from those available for the well-known inverse conduc- 
tivity problem. In fact, on one hand, we prove, Theorem 13. 21 a qualitative form of stability, when a 
very weak topology is assigned on the class of tensors, namely the topology of G-convegence. And, 
on the other hand, we show that a Lipschitz stability bound holds for the mean value of a. 

2 Uniqueness 

For any given K > 1 we consider the class of tensors 

M K = {a G L°°{n,M 2x2 ) | K- 1 ^ 2 < cr£ ■£ < K\£\ 2 for every (el 2 }, (2.1) 
here M 2x2 denotes the set of 2 x 2 symmetric matrices. Let us also introduce 

Sk = W G M K I diver = 0} , (2.2) 
where the null divergence condition div a — is meant in the weak sense 

/ crVu = for every v G H^(Q.) . (2.3) 
Ja 

Let us also introduce 

M = U K >iM K , E = U K >^ K ■ (2.4) 
Given a W 1 - 2 mapping <i> : £1 — > D C M. 2 we denote its Jacobian matrix as follows 

We recall that $ is said to be quasiconformal if, for some Q > 1 it satisfies 

||£>$|| 2 < QdetL»$ a.e. in , (2.6) 

and it is invertible, see for instance Ahlfors pQ. Here, for any matrix A, we denote ||A|| 2 = tr^4^4 T 
and the suffix T denotes transpose. 

For any tensor a G Mk and any quasiconformal mapping $, we introduce 

Tzo-(y) = dct£>$ (.^(y)) , for every y G D . (2.7) 

This new tensor defined in D, called the push-forward of a by is again symmetric and satisfies the 
cllipticity condition with some possibly new K > 1. Moreover, one can verify that such operation 
preserves the bilinear Dirichlet form associated to a, that is 

/ crVu-Vv= / T^ctV^o^- 1 ) • V^o*- 1 ) for every u,v G if 1 (ft) . (2.8) 
Jn J d 

Theorem 2.1. A CT uniquely determines a among all tensors in E. 
The proof will be a consequence of the following two Lemmas. 
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Lemma 2.2. Let a G S. XTiera i/ie tensor 

T^aiD^f (2.9) 

is divergence free in D. 

Proof. Condition (|2.3p can be rewritten as 

/ aVxj • Vw = for every v G ffo(fi) i J = 1, 2 . (2.10) 

By (I2.8[) one computes 

/ T$ctV(($" 1 ) j ) • V(v o = for every v G Hq(Q) , j = 1,2 , (2.11) 

and the thesis follows. □ 

Lemma 2.3. Let a G M and suppose that for a given mapping $ we Ziawe divT$er = in D. 
Then 

div(aL>$ T ) = . (2.12) 

Proof. The proof follows immediately from Lemma l2.2i just by reversing the roles of $ , <& -1 and 
of (j, T$cr, respectively. □ 

We are now in a position to prove our main result. 

Proof of Theorem \2.1\ By the results of Astala, Palvarinta and Lassas [H Theorem 1] , which have 
extended to the L°° setting those of Sylvester [TB] and Nachman [T3] , we have that A a determines 
uniquely the class 

E a = {a 1 G M\a' — T$a, with <3? : f2 — > fi quasiconformal and such that $\qsi — 1} . (2.13) 

The class E a contains at most one divergence free element. In fact, if $ is a quasiconformal 
mappings which fixes the boundary, and such that divider = then, by Lemma l2.3i we have 



div(aD$ T ) = inQ, 

$ = 7 on dn. { > 

Note that this system, is formed by two uncoupled Dirichlct problems for the two components of the 
mapping $. On the other hand we observe that, if a is divergence free, then the identity mapping 
I is itself a solution to (|2.14[) and by uniqueness for the Dirichlet problem, we obtain <I> = I on fl 
□ 

Remark 2.4. It is worth mentioning, that the crucial fact used in this proof is that all linear 
functions are solutions to the elliptic equation (|1.3|) . this is a condition on a which is in fact 
equivalent to the null divergence condition (I2.3[) . Indeed this property has been used already in 
a study on optimization of tension structures in the different context of variational problems and 
G- convergence, see ffjjj, ^ and also Section^ below. 

A further application of this property of linear functions is the possibility to identify the traction 
T applied on the boundary. For every ( 6 I 2 let <fi^(x) = £ • x be the Dirichlet data. Since the 
corresponding solution is u^{x) = £ • x all over 17, we have 

A CT c/j£ = erVii£ ■ v = ■ u — av ■ ^ — T ■ £ , (2-15) 

whose knowledge for every £ G K™ is equivalent to the knowledge of T . 

If the same argument is applied to the particular case of a square network Q = [0, 1] x [0, 1] , 
that is, a portion of fabrics made by two families of parallel elastic strings which cross orthogonally, 
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the identification of a is immediate. The particular situation leads to define as admissible all the 
tensors of the form 

tri{x 2 ) \ 
02(2:1) J ' 

with K^ 1 < ci,(T2 < K. Consider the Dirichlet data <p(x) = x\ on dQ, so that also u(x) = x\ on 
Q, and the corresponding ip{ x 2) = A ff <p on the edge x\ = 1 where v = (1, 0). Then we have 

ci(x 2 ) = <t(1, x 2 )Vu ■ v = ip(xi) , 

likewise, 02(2:1) can be identified as well. 

3 Stability 

This section is devoted to the continuity properties of the inverse of the map 

when we assign to the topology of G-convergence. Let us recall here the basic notions and 
some important properties of the G-convergence. A wide literature is available on this subject, we 
refer for example to the classical papers [TQl [13l [16j ttZ] an d to the book by Dal Maso 9] where 
G-convegence is cast in the more general theory of T-convegence. 

Definition 3.1. A sequence {ah} C Mr is said to G-converge to a £ Mk, and we write 
ah — > a, if for every f £ H^ 1 ^) the corresponding sequence {uh} C Hq(Q) of solutions to the 
inhomogeneous problems 

— div(CT; l Vw/ l ) = / in fl , Uh — on dfl , (3-1) 

converges weakly in H^fl) to the solution u £ Hg(Cl) of the problem 

— div(crVu) = / in O , U = on dft . (3-2) 

It is well known that G-convergence is induced by a compact metrizable topology on Mk, |17[ 
Remark 4] . 

It is also worth recalling that the Lj^-strong convergence implies the G-convergence, [THl Propo- 
sition 5], [17, Remark 11]. 

Theorem 3.2. Given K > 1, the mapping Ejf 3 a — ► A CT £ if (H 1 / 2 ^), F _1 / 2 (0)) has a 
continuous inverse when T*k is endowed with the topology of G-convergence. 

Proof. First we recall that [7] , by the characterization of Y^k as the subclass of those a £ Mk for 
which all linear functions are solutions to (|1.6[) . implies that is a closed set in the G-topology 
and hence it is compact. Let {ah} C £k and a £ be such that ||A CTjl — A CT || — > 0. By the above 

mentioned compactness , there exists a subsequence {a Th } of {ah} such that a Th — > <r' £ and 
we prove cr' = cr. For any ip £ 7J 1 / 2 (9f2), let Ufe, u' the solutions to (|1.6[) when a is replaced with 
with ct/j, cr', respectively. Then, by convergence of the energies [IT] . we have 

< A CTr ip, ip >— / a rh \7u rh ■ Vu rji — > / cr'Vu' ■ Vu' =< A a /ip, cp > , (3-3) 
Jn * Jn 

on the other hand, A ah — > A CT , therefore < A a /tp, tp >=< A CT <y9, </? > for every 99 £ i? 1 / 2 (9J7). From 
the uniqueness Theorem 12.11 we get cr' = a. The above argument applies to any subsequence of 
{ah}- Thus we have obtained that for any subsequence of {ah} there is a sub-subsequence which 
G-converges to cr, and hence, the full sequence {ah} must G-converge to a. □ 
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Remark 3.3. It is worthwhile to compare this result with the case of the inverse conductivity 
problem, that is when the unknown a £ Mk is a-priori known to be isotropic, that is a = 7/ where 
7 G L°°(fi) is a scalar function satisfying K~ x < 7 < K and I denotes the identity matrix. Indeed 
for the inverse conductivity problem, stability with respect to G- convergence fails, see JlSf for related 
arguments. In fact, as is well-known, Marino and Spagnolo J13f proved that there exist a constant 
c > 1, depending only on the space dimension n (in our case n — 2), such that any tensor in M K / C 
can be approximated in the sense of G-convegence by isotropic tensors of Mk- Hence if we had 
stability with respect to G-convergence for isotropic tensors, that would imply the uniqueness in the 
class Mx/c °f anisotropic tensors, which, by the above mentioned example of Tartar cannot hold 
true. Hence the stability result above, Theorem \S.S[ is crucially based on the property of our set of 
admissible matrices, T,k, of being G- closed. 

We recall that in [BJ, and in [7] in a more general context, it was proved that on E^- the G- 
convergence is equivalent to the L°°(r2)-weak* convergence. Therefore, as a consequence of Theorem 
13.21 we also obtain that for every ip £ and for every i,j = 1, 2 the functional F defined by 

F(a) = [ ^o-ij ,<tgS k , (3.4) 
Jn 

depends continuously on A CT . 

In the very special case when, in f|3.4[) . we choose ip = t^t a concrete stability estimate can be 
obtained. In fact, in the next Proposition we show that the average of a, a quantity which can be 
interpreted as a global measure of the pretraction field, depends in a Lipschitz continuous fashion 
on the Dirichlet-to-Neumann map. 

Proposition 3.4. For any a, a' € £ we have 

< (1 + (diam n) 2 ) || A a - A(j/ 1| . (3.5) 

Proof. Being linear functions solutions, we can use (|1.9p . with u = Xi , v = Xj ,i,j = l, 2, both for 
a and a'. We obtain 

(a - a')ij =< (A a - Av>)xi, Xj > , (3.6) 

Jn 

and the thesis follows by straightforward computations. □ 

Remark 3.5. Also in this case it may be interesting to make a comparison with the inverse 
conductivity problem. In fact it is an open problem whether, for the average jij f n 7 of an isotropic 
tensor a = jl, the Lipschitz stability in terms of the corresponding the Dirichlet-to-Neumann map 
holds true, see J3jj. 
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